Neutrino and antineutrino oscillations from decay of Z° boson 
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The two-particle wave function of neutrino and antineutrino after the decay of a Z° boson is 
found as a solution of an initial value problem. The wave packet of this state involves the coherent 
superposition of massive neutrinos and antineutrinos resulting in independent oscillations of their 
flavor counterparts. Possible experiments to observe oscillations of a single lepton or both leptons 
from the same decay or the independently emitted pairs v^v a , v^v^ and v T v T are considered. 
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The standard theories of massive neutrino oscillations 
consider the time-evolution of a neutrino independently 
from a recoil. In spite of the fact that almost all of these 
phenomenological theories provide the canonical formula 
for the probability of oscillations, some important ques- 
tions arise which are discussed in Ref. It has been 
shown in my recent work [2] that, in a two-particle decay, 
noninteracting particles, a neutrino and a recoil, do not 
evolve separately due to quantum correlation between 
the space-like separated objects. Because of this the de- 
cay energy converts to the kinetic energy of the recoil- 
neutrino pair. Thus, I have rigorously proved that the 
neutrino mass eigenstates produced in the electron cap- 
ture decay have the same energy. The question remains: 
does this result still stand for other two-body weak de- 
cays? 

In this communication I apply the formalism developed 
in Ref. Q to the description of the decay of a Z° boson, 
Z° — > v a +v a , where v a is the neutrino with a fixed flavor 
(a = e, /i, r). In spite of the exotic source and the long os- 
cillation length compared with the Earth-Moon distance 
this problem has physical interest because two entan- 
gled particles, neutrino and antineutrino, may oscillate. 
This is yet another peculiar realization of the Einstein- 
Podolsky- Rosen thought experiment 3]. The question 
as to whether neutrino oscillations are observable in the 
invisible Z°-decay has been considered by Bilenky and 
Pontecorvo [3] and Smirnov and Zatsepin However, 
the consensus among these authors on this subject does 
not exist. Let us note also that Z°-decay differs notice- 
ably from the electron capture decay because of the small 
parameter a — Q/M (where Q is the decay energy and 
M is the parent particle mass) used in Ref. [2[ is equal 
to 1. 

The analysis presented below is based on the solution 
of the time-dependent Schrodinger equation with initial 
wave function of a decaying Z° boson taken in the form of 
a finite-size wave packet in its rest frame. The Weisskopf- 
Wigner theory of spontaneous emission @ is used to 
solve the coupled system of differential equations describ- 



ing the time evolution of the boson and the neutrino- 
antineutrino pair in the course of the decay process. The 
wave function which describes the spatiotemporal behav- 
ior of the pair is obtained at times t > 1/r, where T is 
one half of the rate of the decay. In the far zone approx- 
imations this function is represented as a product of the 
relative and the center of mass wave functions, but it does 
not factorize in the neutrino, r„, and antineutrino, r p , 
spatial coordinates. Thus neutrino and antineutrino are 
entangled. It will be shown that the coherent superpo- 
sition of the neutrino and antineutrino mass eigenstates 
has a fixed kinetic energy equal to the decay energy. 
The decaying system is described by the Hamiltonian 



(1) 



where p is the boson momentum operator, Hq is the un- 
perturbed Hamiltonian of the boson, neutrino and an- 
tineutrino fields, respectively, and the weak interaction 
Hamiltonian is given by 



M fa (k,p) 



u: 



iU a jb pC ip / 2+k dj p/2- 



k+h-c, (2) 



where U is the mixing matrix, b, c and d are the second 
quantized operators of the boson, the three (i, j = 1,2,3) 
massive neutrino and antineutrino, respectively. The 
transition matrix element from the initial state a to the 
final one / has the form 



Mfa- 



--g z e a u L (k+p/2)Y 



; (l-7 5 )^(-k+p/2), (3) 



where gz is the coupling constant, Ul and vr are the 
Dirac spinors for the left-handed neutrino and the right- 
handed antineutrino, and the Z°-boson transverse po- 
larization e a = — (0, 0)/v2 corresponds to the state 
in which its spin is directed along the axis z of the rest 
frame. 

The time-dependent perturbation theory is used to de- 
termine the temporal evolution of a decaying state. The 
solution of the Schrodinger equation with the Hamilto- 
nian (JTJ) is sought by using the following ansatz for the 
wave function 
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5^B iJ (p,k,t)ct p/a+k d J - p/2 _ k |o)t 



-iEijt 



(4) 



where |0) is the vacuum state and the eigenvalues of the 
non-perturbed Hamiltonian are 



Ea = M+ 2Fr Eij = e «( k + p/ 2 )+^( k -p/ 2 ) ; ( 5 ) 



where £j(k) = \Jk 2 + m 2 is the energy of the massive 
neutrino or antineutrino with the mass TOj. The coeffi- 
cients Bij(p, k, t) can be found in Weisskopf-Wigner ap- 
proximation just as they where in Ref. Q 



B i:j (p,k,t)-- 



A Q (p)M fa U* l U a 



i(Eij~E a )t-Tt 



>2M{E t3 -E a + iT) 
where the initial value of the coefficient A, 

A (p) = (2^d) 3 / 2 exp(-~d 2 p 2 



(6) 



(7) 



is the boson wave function in the momentum representa- 
tion which describes boson spatial localization before de- 
cay. The one half of the rate of the transition Z° — > v a +P a 
is given by 

r =ibY,\ U ^ 2 l^l 2 \M fa {k,p)\ 2 5(E t] -E a ), (8) 



2M 

The small dependence of T on the boson velocity is ig- 
nored because p <C k, and ultrarelativistic approximation 
reduces this expression to a conventional value 



9%M 

1927T 



86 MeV. 



(9) 



Now one can write the spatial part of the neutrino- 
antineutrino wave function in the coordinate representa- 
tion. With Bij taken from Eq. ([6]) at times t > 1/T, it 
has the form 



1 



>2M 



E 



A (p)M fa (k,p)U* t U aj \^) 
e;(k+p/2) + €j (k-p/2)-M-^ + iT 



x exp{ikr+ipR-i[e i (k+p/2)+e i (k-p/2)]i}, (10) 

where r = r„ — r p is the coordinate of relative motion 
(RM) and R = (r„ + r p )/2 is the center of mass (CM) 
one. 

This function is inconvenient to consider the entangle- 
ment and neutrino oscillations, because it involves inte- 
gration over k and p. The analytical integration over 
k can be performed only approximately. The key ap- 
proximation is referred to as the far zone approximation, 
which means that kr ~ Mt > M/T ^> f. In this ap- 
proximation the main contribution to the integral over 
dfik is given to those k close to the direction of r. In 
line with this assumption, one can put k parallel r every- 
where in the integrand of Eq. (fTU|) except in the factor 



exp(zkr). The remaining function is easily integrated to 
give two terms proportional to exp(zfcr) and exp(— ikr) 
corresponding to outgoing and incoming spherical waves. 
In the far zone, the incoming wave gives an exponentially 
small contribution which can be neglected Q. 

For the outgoing wave, the vectors k and r are paral- 
lel and their direction is determined by the unit vector 
n(0,<p) = r/r, which is the characteristic feature of a 
far zone. Thus stated, the result of the integration of 
Eq. (TlT)|) over dfl^ is given by 



tf(R,r,t): 



(27^) 2 ^/2Mr 



J2Mp)e lpR U* m U aj \is t D 3 } 



expjifcr — i (fc)+p 2 /4e i:) (fc)]t}fcdfc 
2e 4i (k) — M +p 2 /4e 4J (k) - p 2 /2M + iT ' 



(11) 



where the small parameter p/M has been used. The ma- 
trix element (|12p in this approximation transforms to the 
form 



M fa (n) 



gzM 
2 v / 2 



(I — cos ( 



(12) 



The energy t%j{k) depends on the mass term (mf +to 2 )/2. 

The remaining integral in (ITTI) is calculated using the 
residue method by changing the variable of integration 
from k to e and expanding the former around the half of 
the decay energy M/2 



k(e) = h 



-(e-M/2), 



(13) 



where fcjj and the group velocity of neutrinos, — 
2kij/M, are taken at the energy M/2. The lower limit 
of integration in (TlTT) can be extended to -co (since the 
contribution to the integral falls of sharply with increas- 
ing |e| owing to M > T) and the integrant is continued 
analytically in complex plane e. In the lowest order of 
the small parameter p/M the pole of the integrand is 



e p = (M - »r)/2, 
and the integral is evaluated straightforwardly 

k ij (e)e iP ^de _ 2mh 



(14) 



Vij(e)(e-ep) 



^e lPo O(2v l3 t-r). (15) 



Here the contour C encloses the pole (fTl)) . The unit 
step function is due to the requirement of the integral 
convergence on the contour C. The exponent P(e) is 
defined as follows 

P(e) = [kij + (e- Af/2)/« tf ] r- (2e + p 2 /4e) t, (16) 

and its value in the pole can be transformed to the form 



Pn 



-Mt + kar 



p 2 t iT 



2vi.it— r) 



(17) 
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Now, the wave function of relative motion can be con- 
structed by using Eqs. (fTTj) , (fT5"j) and (ITT)) . In doing so, 
it is convenient to express the coupling constant gz in 
terms of T @ . As a result the RM wave function can be 
written as 



where tcM = Md 2 is its spreading time. This time is less 
than the flight time if d < 10~ 3 cm. 

The total wave function describing the evolution of the 
neutrino-antineutrino pair after decay is a product of the 
CM and RM parts 



^(n,r,t) = F(n)Y,R ij {r,t)U* i U aj \v i D j ), (18) $(R,r,f) = f CM (R,^(n,r,<)expHMt). (25) 



'■j 



where its radial part involves the functions 

2vid — r 



Vr 

r 



Rij (r, t) =^f- exp ( ikijr - " ' )Q(2v ij t-r), (19) 



depending on r and t, and the angular one, 



cos( 



(20) 



depends on the orientation of the vector r. The func- 
tion (|18p is normalized in the ultrarelativistic limit, when 
R oe ik ii r where 



Vij = 1 and Rij 



exp [-T(t-r/2)} Q(2t- 



(21) 



It is seen that the energy of this state is equal to the 
decay energy, Q = M. This is a general feature of the bi- 
partite decay into noninteracting fragments constrained 
only by momentum and energy conservation. The wave 
packet | ^>(R, r : t) | 2 increases with time in a longitudinal 
direction (along the vector r) with the velocity 2 and in 
a transverse one with the velocity l/(27r 3 / 2 Md). Hence, 
the spatiotemporal behavior of the joint quantum state 
of the neutrino and the antineutrino following Z°-decay is 
in agreement with the results obtained in the theoretical 
studies of the electron capture decay [2j . The distinctive 
feature of our system is the coherent superposition of the 
two leptons with flavor mixing, which both oscillate. To 
describe this phenomenon, it is necessary to consider the 
two-particle wave function (|2"5")) in the observable coordi- 
nates of the neutrino and the antineutrino. This normal- 
ized function has the form 



The function tp describes a coherent superposition of ex- 
ponential wave packets of massive neutrino-antineutrino 
pairs with different sharp edges at r = 2vijt and widths 
Dij = 2vij/T ~ 10~ 13 cm. The later depend only on 
the dynamics of the decay process and may change later 
due to dispersion. To estimate such spreading effect, it is 
necessary to extend the series expansion of the function 
k(e) in Eq. (fT3")) up to the second order in e — M/2. This 
gives rise to an additional integral in Eq. (fTTj) as shown 
in Ref. [2J . The calculations described in this work gives 
the spreading time of the RM wave packet 



tRM 



M 3 Df. 



8vij(m 2 + m|) 



100 s. 



(22) 



The time required for a neutrino or an antineutrino 
to propagate the distance comparable with oscillation 
length, "flight time", is 10 times smaller than t^M- This 
allows to use the sharp edge wave packets for analysis of 
neutrino oscillations in Z°-decay. 

The integration over dp in Eq. (jlll) can be readily per- 
formed by using expressions ([Jj and (fT7|) . The resulting 
CM wave function has the standard form for the decaying 
bipartite system 



*ca/(R,0 : 



1 



r 3/4 



id' 



_*t\3/2 
Md) 



exp 



R 



(23) 



This is a spreading wave packet with the time-dependent 
width 



D R (t) = 



1 V2 



d, t <c ten 
t/Md, t > t C M 



(24) 



CM 



t )i)(n,Y v -Y i) ,t) e 



-iMt 



(26) 



The function does not factorize in these variables - a 
direct indication of the spatial entanglement of leptons. 
Each of the three massive neutrinos is entangled with 
three massive antineutrinos. 

It is instructive to consider the recoilless decay of an 
infinitely heavy mother particle (M — > oo). Setting p = 
in Eq. ([S]) and performing exactly the same calculations 
as above with A$ — 1, one can obtain the function of the 
pair 



iML 



(27) 



This state is also entangled because massive neutrinos 
and antineutrinos are not emitted in a momentum eigen- 
state due to the fact that the interaction Hamiltonian @ 
commutes only with the total momentum. 

The spatial part of the neutrino-antineutrino wave 
packet (|26l) for fixed t is proportional to a product of 
the Gaussian and exponential functions 



,!, m2 I (r„+rp) 2 |r„-r p | 
\v(T v ,ri?)\ ~exp<i 7-^ h2- 



D 



(28) 



where Dij ss D = 2/T in the ultrarelativistic approxima- 
tion. The probability distribution of the two-particle co- 
ordinates depends on r„, r p and the angle between these 
vectors. Therefore, the packet has an axially symmetri- 
cal shape with respect to the axis passing through the 
CM in the direction of the vector n. It is easily to see 
that probability reaches its peak when the vectors have 
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equal lengths and are located on the symmetry axis an- 
tiparallel to each other. Hence, the daughter particles 
are traveling back-to-back in the rest frame of Z° as pro- 
vided by the classical picture. However, the decay is not 
isotropic because of the spin polarization of the boson. 
The preferred direction in which the particles travel is 
determined by the factor ^(n)! 2 ~ (1 — cos6*) 2 , which is 
to say the antineutrino travels mainly in the direction of 
the boson spin. 

The structure of the neutrino- antineutrino wave packet 
of the function (|26|) suggests two kinds of experiments to 
observe flavor oscillations. 

Neutrino and antineutrino oscillations in one detec- 
tor experimen. The first is a well-established experiment 
for detecting neutrino or antineutrino oscillation. Sup- 
pose that only a neutrino is detected regardless of the 
antineutrino position. In the detection of a flavor neu- 
trino, when the position of a detector is scanned, such 
measurements allows to determine the oscillation picture. 
The flavor-changing transition v a — > vp is determined by 
the probability density 



dr v 



dr.- 



i 



(29) 



Changing the variable of integration from r p to r, one 
gets the integral 



dr 

— exp 



(r+2r, y 



2r 
~D 



where 



\F{6^)\ 2 W a p{r)Q{2t-r), 
(30) 



W aP {r) = S a p - 4j2U ai U m U aj U Pj sm 2 (<^(r)) (31) 

is the probability of the neutrino flavor transition v a — > 
Vp in the case of the CP invariance and 



<t>ij (r) 



i>>~ 



2M 



(32) 



is the oscillation phase. The integrand of Eq. (|30|) has a 
sharp maximum at r — 2t, and integration can be done 
approximately since the width D of the exponential func- 
tion is much smaller than that of Dn(t) of the Gaussian 
one, if Ft ^> 2d/Xz, where Xz is the Compton wave- 



length of the boson (r ~ 10 



23 



The integration is 



straightforward, albeit somewhat tedious, and its result 
is 



dP a 







dr,,dVl, 



= W a p(2t)(l - cos d v + cos 2 ^) 



3r 2 



A^/ 2 D\{t) 



exp 



(33) 



where is the angle between the outgoing neutrino and 
the direction of the Z° boson spin. The probability den- 
sity of the oscillation transition D a — ¥ Dp may be derived 



in perfect analogy to above calculations, and has the form 



dP a p 
dr^d^ljj 



■ W a p{2t){\ + costfp + cos 2 tf p ) 



3r 2 



4ir 3 / 2 D 3 R {t) 



exp 



lit) 



(34) 



because W a p — W a p. 

It is seen that the probability densities (1331) and (fM)) 
involve the amplitude of oscillations depending on the 
spatio-temporal localization of a neutrino. It owes its 
origin to the evolution of an unstable localized system 
during its decay. Obviously, the amplitude can not be 
found without rigorous treatment of the decay process as 
described above. The probability of transition is obtained 
by integration of (j3"3")) or ([54]) over a spatial coordinate, 
and it is equal to W a p(2t). 

The detection time t is equal to distance between the 
detector and the source, t = r v = r/2. Therefore, the 
oscillation length, associated with Am- ■ , is 



= 2tt- 



M 



(35) 



where the neutrino energy is obviously equal to M/2. 
Two conditions are necessary for the neutrino oscilla- 
tions to be observed. The first condition is the coher- 
ence of different mass eigenstates. Coherence is pre- 
served over distances not exceeding the coherence length 
L C oh- The latter is defined as the distance at which the 
phase difference due to energy spreading of a neutrino, 
Se = r, obeys the equation (dcj)/de) m 5e = 2ir, which 
gives L co /j = L osc Q/T ~ 10 L osc . The second condition 
requires that at a distance, which is comparable with 
the oscillation length, the spatial separation of the wave 
packets corresponding to different mass eigenstates would 
be considerably less than the packet width. Such, indeed, 
is the case 



L osc Sv TL C 



vD 



2v 2 



d 



dm 2 ) 



Am 2 = 7T— < 1. (36) 



Hence, the spatial separation of the massive neutrino 
wave packets does not affect the observation of neutrino 
oscillations in any two-particle decay. 

All of the above applies to the decay of a single Z° bo- 
son. In the case of a steady flux of Z 0, s whose decays pro- 
vide the incoherent fluxes I a of the three types of neutri- 
nos (a = e, /i, t), the intensity of neutrino /3 at an any dis- 
tance from source is proportional to ^ Q I a W a p(r). Be- 
cause the probabilities of emission of different neutrino- 
antineutrino pairs are equal, I a — la, the sum is equal to 
Io according to Eq. (|3"Tj) . Hence, the neutrinos (as well 
as antineutrinos) do not oscillate, which is in agrement 
with the result of Ref. Q. 

Neutrino and antineutrino oscillations in coincidence 
measurements. The term "coincidence" means that the 
flavors of both leptons from the same decay are fixed by 
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two detectors. The flavor-changing process v a v a 
is determined by the probability density 



dP - 

Uj± aa- 



dv^dvp 



(37) 



By using Eq. (|2"Tj) and (|3ip . one finds 



*cm(- 



dYudYu 



x \F(eip)\ 2 R 2 (r,t)W af} (r)W aj (r) 



(38) 



The highest output is obtained when the detectors are 
located on the symmetry axis at equal distances from the 
source r„ = r p = t. For this configuration, the oscillation 
amplitude is 



(1 -cosOfTrl 
647r 5 /2 j D|( rv ) 



(l-cos6Q 2 r (d_ 
64ir 5 / 2 r„ \X Z 



(39) 



This amplitude corresponds to the highest spatial corre- 
lation of the two leptons at a distance equal to 2t, that 
is in complete agreement with the classical picture of a 
bipartite decay. 

The probability of detecting the neutrino vp and the 
antineutrino v 1 from the decay Z° — > v a + v a is obtained 
by the integration of Eq. 



(40) 



It is seen that a neutrino and an antineutrino oscillate in- 
dependently, and their oscillation lengths are determined 
by Eq. Q36p. If the pairs v e v ei v^v^ and v r v T arc emitted 
independently, and their fluxes are equal, then the prob- 
ability of observing z/g in one detector and T> 1 in another 
one will be equal 



Was^pf = dP a& ^ = W aP {2t)W ai {2t). 



W M =^W af) W ai , 



(41) 



due to mixing of signals from two detectors. In order 
to observe the independent oscillations of neutrinos and 
antincutrinos, it is necessary to synchronize the time on 
two detectors to make sure that the leptons come from 
one decay. 



In summary, an accurate analytical solution for the 
joint quantum state of flavor neutrino and antineutrino 
following the decay of a Z° boson has been found. The 
evolution of the state provides yet another exactly calcu- 
lable illustration of the famous Einstein-Podolsky-Rosen 
thought experiment. The new effect is the entanglement 
between two oscillating leptons. It is shown that the 
oscillations of neutrino and antineutrino from the same 
decay proceed independently and may be observed by 
measuring the flavor of both particles. Massive neutrinos 
and antincutrinos are not emitted in a state with definite 
momentum, however the energy of the two-lepton wave 
packet is fixed and is equal to the decay energy. Hence 
it can be argued that the neutrino mass eigenstates pro- 
duced in bipartite decay have the same energy. 



Despite its purely academic interest, the solved prob- 
lem is essential to the understanding of the physics of 
neutrino oscillations. Besides, the sources of neutrinos 
from the Z°-decay, as well as detectors recording neutri- 
nos from a single decay may well be conceivable in the 
future experiments. Accordingly, I have considered the 
experiments allowing to observe oscillations of one lep- 
ton as well as both leptons from the decay of a single 
boson. If three neutrino-antineutrino pairs are emitted 
independently, the oscillations can be observed also by 
the two detectors connected into the coincidence circuit 
with one output or synchronized in time. A similar result 
was obtained in Ref. [a]. However, the results of that 
study are based on the erroneous assumption that Z° 
boson decays into a coherent superposition of neutrino- 
antineutrino pairs with different flavors Q. 



which is different from Eq. (1401) for the decay of a single 
boson. For the detectors connected to a coincidence cir- 
cuit with one output, the oscillation pattern is distorted 
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